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A spin model that displays inverse melting and inverse glass transition is presented and analyzed. 
Strong degeneracy of the interacting states of an individual spin leads to entropic preference of 
the "ferromagnetic" phase, while lower energy associated with the non-interacting states yields a 
"paramagnetic" phase as temperature decreases. An infinite range model is solved analytically for 
constant paramagnetic exchange interaction, while for its random exchange, analogous results based 
on the replica symmetric solution are presented. The qualitative features of this model are shown 
to resemble a large class of inverse melting phenomena. First and second order transition regimes 
are identified. 
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We all tend to associate order parameter with order, 
namely, with less entropic microscopic realizations. This 
is indeed the general situation in nature: crystals are 
more ordered than liquids, ferromagnets have less en- 
tropy than paramagnets. Even the entropy associated 
with a glass, an out of equilibrium, frozen frustrated 
state, is less than that of a liquid phase of the same ma- 
terial. 

There are, however, exceptions, where an "order pa- 
rameter" does not reflect order, and the entropy growth 
during crystallization or freezing. The prototype of these 
phenomena is inverse melting, i.e., a reversible transi- 
tion between a liquid phase at low temperatures to a 
high temperature crystalline phase, observed in He and 
He 4 at extreme conditions (temperature below 1°K, pres- 
sure above 25 bar) []]. A similar phenomenon was ob- 
served recently at room temperature and atmospheric 
pressure in P4MP1 polymer solutions Q- Ferroelectric- 
ity in Rochelle salt is another example, where the spon- 
taneous polarization is lost below the (lower) Curie tem- 
perature, this time the transition is second order in type 
[3- The pinned-crystalline inverse transition of vortex 
lines in the presence of point disorder at high tempera- 
ture superconductors 4] is also considered as an example 
of inverse melting. However, in that system, the inten- 
sive order parameter (bulk magnetization) is lower in the 
crystalline phase, and the response functions are higher, 
i.e., the disordered phase is stiffer than the ordered phase. 

Even if the crystalline state is the thermodynamically 
preferred one, the dynamics of the system may prevent its 
appearance. In glass forming materials ergodicity break- 
ing takes place at a finite temperature and the system is 
trapped into a frozen disordered state. One expects that 
an "inverse" glass transition phenomenon, analogous to 
inverse melting, may also take place. An interesting ex- 
ample in polymeric systems is the reversible thermoge- 
lation of Methyl Cellulose solution in water When 
a (soft and transparent) solution of Methyl Cellulose is 
heated (above 50°C, for a 10 gr/liter solution) it turns 
into a white, turbid and mechanically strong gel. Unlike 
the boiling of an egg that involves an irreversible transi- 
tion from a metastable to a stable state, this transition is 
reversible upon cooling, and the polymer is redissolved on 



subsequent cooling. In its high temperature phase, the 
Methyl Cellulose gel exhibits, like many other gels [||, 
glassy features. Non monotonic temperature dependence 
of the glassy order parameter has been already reported 
for a random heteropolymer in a disordered medium ; 
this may be considered as the glassy analogue of the flux 
line crystallization £| . The liquid- liquid transition theory 
for polyamorphous materials predicts an inverse freezing 
transition even for the most known liquid, water. In the 
hypothesized phase diagram presented in |8( a low den- 
sity liquid (at about 150 bar, — 100°C) becomes a low 
density amorphous ice upon heating. 

In many branches of statistical physics the presenta- 
tion of a simple spin model (Ising, Potts, and SK models, 
for example) turns out to be a very beneficial step that 
yields both physical insight and quantitative predictions. 
In this paper, such a model for inverse melting is pre- 
sented and analyzed for homogenous and heterogenous 
systems in the mean field level. The model exhibits both 
inverse melting and inverse glass transition, and allows 
first order and second order transitions. We believe that 
this generic model is applicable for the qualitative de- 
scription of the above mentioned phase transitions (ex- 
cept for the inverse melting in superconductors which 
requires a different model). 

Let us begin with an intuitive argument focusing on 
one of the above mentioned systems, namely, a single 
Methyl Cellulose polymer chain in water. In order to 
explain the inverse freezing it seems plausible to assume 
that its folded conformation is favored energetically while 
its unfolded conformation is favored entropicaly [See fig- 
ure ■ The entropy growth of the open conformation 
may be related to the number of possible microscopic con- 
figurations of the polymer itself, but it may be attributed 
also to the spatial arrangement of the water molecules in 
its vicinity |9(. 

The main cause for inverse freezing is that the " open" 
conformations of the polymer are also the interacting 
structures, as they allow for the formation of hydrophobic 
links with other polymers in the solution, a process that 
leads to gelation. This seems to be a general prescription 
to both inverse melting and inverse glass transitions: the 
noninteracting state is favored energetically, while the in- 
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FIG. 1: A sketch of the energy and entropy dependence on the 
linear size of a Methyl Cellulose polymer in water. The folded 
conformation costs less energy due to more favorable interac- 
tions between hydrophobic sequences along a single chain, but 
are less entropic as water molecules has to arrange in cage like 
structures around the hydrophobic constituents of the chain. 
The unfolded conformation admits much more microscopic 
configurations. The interaction with other polymers in the 
solution is suppressed in the folded state. 
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FIG. 2: Phase diagram and spinodal lines for the ordered 
model Eq. |JSJ in the D — T plane for r = 1 (Blume-Capel 
model, inset) and for r = 6. The value of r — 6 has been 
chosen in order for the effect to be more pronounced, but 
inverse melting is seen for r lower than 2 



teracting state is favored by the entropy. 

Let us now present a very simple model that incor- 
porates these features. It is based on the Blume-Capel 
model 0]>E3 for a spin one particle with "lattice field" 
that lower the energy of the "zero" (noninteracting) 
state. In contrast with the original Blume-Capel model, 
we consider the ±1 spin states (that interact with other 
spins) to be more degenerate. The system consists of a 
lattice of N sites and the Hamiltonian is given by 

N 

H=-jJ2S i S j + Dj2sf (1) 

<i,j> i=l 

where the spin variables are allowed to assume the values 
Si = 0,±1. The summation over < i,j > is over each 
distinct pair once. Turning back to our polymer anal- 
ogy, spin represents schematically the compact non- 
interacting polymer coil, the stretched polymer (interact- 
ing with its neighbors) is represented by spin ±1. The 
positive constant D measures the energy preference of the 
compact spatial configurations, and the "ferromagnetic" 
interaction between spins, J, is related to the concentra- 
tion of polymers (or the pressure). The spin state is 
assumed to be n-fold degenerate, and the ±1 states are 
m-fold degenerate so that r — m/n > 1 is the degen- 
eracy ratio that dictates the entropic advantage of the 
interacting states. It turns out that all the results pre- 
sented here are independent of the absolute degeneracies 
m and n, and depend only on their ratio r. 

Using standard gaussian integral techniques one finds 
an expression for the free energy per spin in the infinite 
range limit: 

/ = F/N = f3Jm 2 /2 - ln[l + 2 r cosh{(3Jm)e~ pD ] (2) 



where m is the order parameter of the system (magneti- 
zation per spin), m = {-k Y2i=i ^i)- The phase transition 
curves are obtained numerically by solving for the mini- 
mum of Eq. J5J) with respect to m. Scaling the temper- 
ature and D with the interaction strength J, the phase 
diagram is shown in Figure I n the inset, results are 
presented for the original Blume-Capel model (i.e., the 
r = 1 case): the line AB is a second order transition line, 
above it is a paramagnetic (m = 0) phase and below it 
the system is Ferromagnetic (m^O). Below the tricriti- 
cal point (B) the phase transition is first order, and the 
three lines plotted are: the spinodal line of the ferromag- 
netic phase BE (above this line the m/0 solution ceases 
to exist), the spinodal line of the paramagnetic phase 
BC (below this line there is no m = minimum of the 
free energy) and the first order transition line BD. Along 
BD the free energy of the paramagnetic phase is equal 
to that of the ferromagnetic state. Clearly, the Blume- 
Capel model displays no inverse melting: an increase of 
the temperature induces smaller order parameter. 

The situation is different as r increases, as emphasized 
by the main part of Figure @ ■ The same phase diagram 
is presented, but now r = 6, so the interacting states 
have larger entropy. The tricritical point is shifted to the 
left, leaving a region of second order inverse melting, and 
the orientation of the BD line also changes, establishing 
the possibility of first order inverse melting. Note that 
the r = 6 transition lines converge to the r = 1 lines as 
T — ► 0, since the entropy has no effect on the free energy 
at that limit. The ferromagnetic phase also covers larger 
area of the phase diagram for r = 6, a fact that reflects 
again its entropic advantage. 

To allow qualitative comparison of our cartoon model 
with experimental results, the appropriate parameters 
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FIG. 3: Phase diagram and the spinodal lines for the ordered 
Blume-Capel model in the interaction-temperature plane with 
r = 6. The interaction J/D represents the concentration 
("pressure"). 

should be identified. There are three parameters in the 
model as it stands: D represents the energetic advan- 
tage of the noninteracting state, r (if larger than 1) is 
the entropic gain of the interacting state, and J is the 
strength of the interaction. In most of the physical sys- 
tems that display inverse melting the controlled exter- 
nal parameter is the strength of the interaction: pressure 
(for He 3 and He 4 ) or concentration of the interacting ob- 
jects (for polymeric systems and Rochelle salt - Ammo- 
nium Rochelle salt mixtures). As long as the only effect 
of the pressure is to increase the strength of the effec- 
tive interaction among constituents, it may be modelled 
by changing J. The resulting phase diagram should be 
compared, though, with the T — J plot of our model pre- 
sented in Figure © . The decrease of the transition tem- 
perature with the increase of interaction strength (pres- 
sure) is physically intuitive, as larger interaction favors 
energetically the ferromagnetic phase. As emphasized re- 
cently by [lj, the slope of the first order transition line 
in the pressure-temperature plane is required by the cor- 
responding Clausius-Clapeyron equation: 

dP _ 02 — Si 

where V%, V\ are the volume (the extensive parame- 
ter conjugate to the pressure) of the solid and liquid 
phases respectively, and S2, Si arc their entropies. In- 
verse melting is possible if the numerator of J3J) is neg- 
ative, so for "normal" transitions (V2 > V\) one ex- 
pects a negative slope of the transition line. In real 
magnetic or electric system the intensive-extensive pairs 
[magnetization-magnetic field (M • dH) or polarization- 
electric field (P • f£E)], appear in the free energy function 
with inverse sign relative to PdV. If the order parameter 
vanishes, or takes smaller values, in the "liquid" (disor- 
dered) phase, this implies also negative slope of the first 



order transition line in the temperature-external field 
plane. An interesting exception is the inverse melting 
of vortex liquid in superconductors, where the magneti- 
zation of the crystalline phase is smaller than that of the 
liquid and the transition line slope is actually positive. 

Inverse freezing, the (reversible) appearance of glassy 
features in a system upon raising the temperature, may 
be incorporated in our model by introducing random cou- 
pling Jij, as in the standard spin-glass models 13] . This 
randomness may fit, in particular, to the gelation tran- 
sition of Methyl Cellulose, as it occur only when the hy- 
drophobic sequences are deposited at random along the 
chain. The random-exchange analogue of the Hamilto- 
nian is: 

N 

H=Y, JijSiSj+^DSf (4) 

<i,j> 2=1 

where the exchange interaction between the i and the 
j spin is taken at random from some predetermined 
distribution. Following the paradigmatic Sherrington- 
Kirkpatrick (SK) analysis [l3| of the infinite range spin 
glass, we assume gaussian distribution of the exchange 
term with zero mean: 

!~N~ NJf, 

where is the width of the distribution. The replica 
trick is then implemented to get the free energy at the 
large N limit. 

The case r = 1, namely the random exchange version 
of the Blume-Capel model, was first introduced and dis- 
cussed by Ghatak and Sherrington (GS) ,14] who used 
symmetric replica to obtain the relevant phase diagram. 
The GS solution seems to display inverse freezing even for 
the r = 1 case, but more detailed analysis by da Costa 
et. al. revealed that the glass order-parameter takes 
nonzero values (with a variety of stability features) in the 
area below the GS transition line, and the temperature 
dependence is monotonic. Recently, the full replica sym- 
metry breaking analysis has been implemented for the 
GS model 17], and the results admit no inverse glass 
transition. Here we present a replica symmetric analysis 
of the same hamiltonian where the interacting states are 
highly degenerate, i.e., r > 1. Following 0], we obtain 
the phase transition and the spinodal lines, and the re- 
sults support, again, both first and second order inverse 
glass transition. 

The replica technique relies on the identity ln[Z\ = 
lim n ^o^{Z n — 1), where Z is the partition function of 
the system and Z n is interpreted as the partition func- 
tion of an n-fold replicated system Si — > Si a ,a — l...n. 
The average free energy may be computed using (3f = 
—lim n ^Q 1 ^(Z n — 1). The disorder average is taken for 
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Z n using the Gaussian distribution JSJ and gives: 
>J 2 



Z n = Tr exp 
(3 2 J 2 



N 



a>b i 



2N 



(6) 



where a,b = X...n denotes the replica. Implementing the 
Hubbard-Stratanovitch identity yields the free energy per 
spin: 



a>b a 



where 



L = 2/3 2 J 2 QabSaSt + /3 2 J 2 J2 Q*aS 2 -fiD^Sl 



a>b 



(8) 



q aa and q a b, the diagonal and the off diagonal entries of 
the "order parameter matrix", are given self-consistently 
by the saddle-point condition: 



q a b — (S a Sb) 
Qaa = 



(9) 



where (...) stands for thermal average. In order to solve 
this model it is necessary to make assumptions on the 
order parameter matrix elements q a b, and the simplest 
ansatz, is symmetry with respect to permutations of any 
pair of the replicas: q a b = q, Va 7^ b, q aa = p, Va. Using 
this replica symmetric assumption one obtains 



1 r, , i\ 

-== \ dz exp 



In [1 + 2 r e 1 cosh(/3J^2qz) 



with 



1 = P 2 J 2 (p-q)-(3D 



(10) 



(11) 



Extremizing the free energy with respect to q and p 
one gets by the following coupled equations: 



dz exp{-^-) Ar 2 e 2 ^smh 2 (f3J^/2q~z) 
2^ [1 + 2re~< cosh(/3JV2^z)] 2 



(12) 



dz exp(-^-) 2re 7 cosh(/3 J^/2qz) 
1 + 2re^ cosh(/3 J*j2qz) 



(13) 



The coupled equations (|12fl and (| 1 31) are numerically 
solved (with the possibility of multiple solutions if more 
than one stable state exists), and the location of the first 
order transition line is then determined by comparison 
of the free energy values (plugging q and p into l |10(l ). 
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FIG. 4: Phase diagram and the spinodal lines for the disor- 
dered model in the D-T plane for a constant interaction J for 
r = 6. 



The resulting phase diagram is shown in Fig. Q for the 
case r = 6, and displays all the essential features that 
exist in the ordered model, including a tricritical point 
and spinodal lines. 

To conclude, the basic theoretical insight of Blumc 
and Capel, to have a spin model with low energy non- 
interacting (zero) state, may yield an inverse melting 
transition once the model is enriched with an entropic ad- 
vantage of the interacting phase. It should be emphasized 
that the higher entropy associated with the interaction 
do not unavoidably entail inverse melting; this property 
may be "buried" below energetic and other constraints 
that dominate the system, yet it may change the phase 
diagram predicted by the naive assumption that higher 
energy implies higher entropy. 

The authors wish to acknowledge Prof. Y. Rabin for 
most helpful discussions of the subject. 
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